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Abstract 

In  this  report  we  will  focus  on  minimising  the  cost  due  to  inspection  and  repair 
of  hydraulic  structures.  The  optimisation  is  based  on  symmetric  properties  of 
the  underlying  physical  deterioration  process:  i.e.  the  damages  per  time-unit  are 
exchangeable  and  the  probabilities  of  preventive  repair  and  failure  are  obtained 
by  conditioning  on  the  average  amount  of  deterioration  with  regard  to  a  finite  or 
an  infinite  time-horizon.  By  introducing  a  prior  for  the  average  deterioration  per 
time-unit  we  can  account  for  uncertainty  in  the  decision  problems.  Advantage  i 
of  our  Bayesian  approach  arc  that  we  base  our  mathematical  models  on  an  ob¬ 
servable  quantity,  namely  the  damage,  and,  surprisingly,  that  our  results  are  just 
sums  of  products  which  can  be  easily  evaluated.  Two  examples  from  the  field  of 
hydraulic  engineering  arc  studied:  determining  a  preventive  repair  interval  when 
a  safety  norm  is  given  and  a  cost-optimal  periodic  inspection  rate  if  there  is  a 
possibility  for  a  preventive  repair  during  each  inspection. 
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1  Introduction 


This  report  will  discuss  the  subject  of  preventive  maintenance  of  hydraulic  structures. 
The  problem  emerges  from  the  field  of  structural  engineering,  although  the  same  ap¬ 
proach  can  be  applied  within  the  field  of  mechanical  and  electrical  engineering. 

Preventive  maintenance  of  hydraulic  structures  consists  mainly  of  two  parts:  in¬ 
spection  and  repair.  A  repair  can  either  be  a  preventive  repair  (before  failure)  or  a 
corrective  repair  (after  failure).  Inspection  is  often  needed,  because  it  is  not  possible 
to  observe  continuously  the  physical  condition  of  a  structure,  e.g.  the  height  of  a  dyke 
or  the  strength  of  a  bridge.  For  an  extensive  problem  description  about  maintenance 
of  hydraulic  structures  see  van  Noortwijk  [11]. 

Inspection  and  repair  entail  substantial  costs.  Money  can  be  saved  by  weighing 
inspection  cost  on  one  hand  with  repair  and  failure  cost  on  the  other.  A  high  inspection 
rate  will  result  in  high  inspection  cost  and  low  repair  and  failure  cost.  No  inspection  at 
all  will  surely  lead  to  failure  and  high  cost.  An  optimal  or  near-optimal  maintenance 
policy  lies  somewhere  in  between.  In  practice  the  failure  cost  is  often  much  higher  than 
the  cost  of  preventive  repair,  where  failure  cost  includes  both  cost  of  corrective  repair 
and  other  public  or  private  losses  due  to  failure. 

To  lower  the  expected  sum  of  inspection,  preventive  repair  and  failure  cost  a  wide 
variety  of  mathematical  models  has  been  built.  For  surveys  see  Barlow  &  Proschan  [4], 
Pibrskalla  &  Voelker  [12],  Sherif  &  Smith  [13]  and  Sherif  [14].  In  most  cases  inspection 
and  repair  are  modelled  in  a  probabilistic  way.  Quite  popular  are  maintenance  models 
based  on  lifetime  distributions  assuming  that  the  deterioration  law  of  a  structure  is 
Markovian.  A  problem  in  practice  is  to  get  data  for  either  a  lifetime  distribution  or 
the  transition  probabilities  of  a  Markov  chaiu  (where  a  transition  is  a  change  of  the 
structure’s  physical  condition).  Unfortunately,  only  a  few  applications  can  be  found. 

Instead  of  assuming  a  lifetime  distribution  or  a  Markov  chain,  we  will  follow  a 
different  approach.  As  iu  structural  engineering  we  will  define  a  failure  as  the  event 
at  which  a  structure’s  resistance  drops  below  the  stress.  We  will  focus  on  stochastic 
resistance  and  a  constant  deterministic  stress. 

The  resistance  will  be  modelled  with  damages  appearing  per  time-unit.  The  ad¬ 
vantage  of  this  approach  is  that  damage  is  observable,  while  a  lifetime  or  a  transition 
probability  matrix  is  not.  Furthermore  we  do  not  make  the  usual  assumption  of  mu¬ 
tual  independent  damages,  but  base  our  proposition  on  the  physics  and  the  decision 
problem. 

Our  Bayesian  approach  is  applicable  for  decision  problems,  where  the  joint  proba¬ 
bility  density  of  the  amounts  of  damage  per  time-unit  is  invariant  under  permutation 
and  depends  on  the  average  amount  of  damage,  while  considering  a  finite  or  infinite 
time-horizon.  These  properties  arc  known  as  exchangeability  and  /j  symmetry  respec- 
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tively.  They  will  be  discussed  in  section  2  and  were  first  applied  to  decision  problems 
by  Mendel  [9,  10]. 

An  advantage  of  the  discretisation  in  time-units  and  the  application  of  exchange¬ 
ability  and  k  symmetry  is  that  the  calculation  of  the  probabilities  of  preventive  repair 
and  failure  becomes  tractable.  All  probabilities  can  be  expressed  in  sums  and  prod¬ 
ucts.  Moreover,  the  Bayesian  modelling  is  based  on  the  underlying  physics.  In  general 
a  classical  Bayesian  provides  a  prior  for  every  uncertain  parameter  in  a  probabilistic 
problem  formulation,  without  accounting  for  ways  to  determine  these  priors  or  even  to 
make  sure  whether  the  parameter  is  observable  in  practice. 

Concerning  to  the  inspections  we  will  focus  on  a  periodic  inspection  policy,  although 
Barlow,  Hunter  k  Proschan  [2]  proved  that  a  constant  inspection  rate  is  not  optimal 
when  the  lifetime  distribution  is  not  exponential.  Our  failure  distribution  based  on 
accumulated  damages  will  not  be  exponential.  In  practice  a  periodic  inspection  policy 
is  often  chosen,  because  it  is  easier  to  plan  manpower.  In  the  future  it  is  interesting  to 
look  at  a  resistance  dependent  inspection  policy. 

Two  examples  will  be  studied  in  section  3. 

First,  we  focus  on  determining  a  preventive  repair  interval  (without  needed  or 
possible  inspection),  while  the  failure  probability  should  not  exceed  a  certain  safety 
norm. 

Second,  wo  study  an  inspection  model  which  is  an  extension  of  the  basic  model  of 
Barlow,  Hunter  k  Proschan  [2].  It  includes  periodic  inspection,  while  at  an  inspection 
a  preventive  repair  only  will  be  executed  if  the  resistance  is  below  a  preventive  repair 
level.  If  an  inspection  is  scheduled  too  late  the  structure  can  pass  into  the  failed  state. 
A  large  amount  of  money  is  involved  when  this  happens.  Passing  the  preventive  repair 
level  can  only  be  noted  through  inspection,  whereas  a  failure  will  be  noted  immediately, 
like  an  inundation  of  an  area  protected  with  flood  defences  against  the  sea.  The  time- 
horizon  can  be  Unite  and  the  decision  variable  is  the  inspection  rate. 

All  technical  proofs  can  be  found  in  the  Appendix. 
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Notation 


cp  Cost  of  failure. 

cj  Cost  of  one  inspection. 

cp  Cost  of  preventive  repair. 

C(k)  Expected  sum  of  inspection,  preventive  repair  and  failure  cost. 

Di  Amount  of  damage  in  the  1th  time-unit. 

kA  Length  of  inspection  interval. 

NA  Time-horizon. 

Pnorm(n)  Safety  probability  norm, 

r  Resistance  at  time  zero. 

R  Resistance  or  strength. 

Rn  Resistance  in  time-unit  n. 

R  -  S  <  0  Failure. 

s  Constant  stress  or  failure  level. 

S  Stress,  load  or  action  effect. 

A  Length  of  time-unit. 

9  Average  amount  of  damage  per  time-unit. 

p  Preventive  repair  level. 
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2  Assumptions  and  definitions 

In  the  field  of  mechanical  and  electrical  engineering  one  often  considers  lifetime  dis¬ 
tributions  to  model  the  appearance  of  failure  in  a  system,  e.g.  a  motor  or  switch  is 
working  or  not,  whereas  in  the  area  of  structural  engineering  failure  is  decomposed  to 
comparing  a  structure’s  resistance  or  strength  R  with  its  stress,  load  or  action  effect  S 
(see  e.g.  Ang  &  Tang  [1]  and  [6]).  A  structure  is  said  to  fail  if  its  resistance  is  below 
the  stress,  i.e.  if  the  so-called  performance  function  R  -  S  is  below  zero.  In  general 
both  resistance  and  stress  are  unknown  functions  of  time  and  moreover  not  necessarily 
independent.  In  this  paper  we  will  restrict  ourselves  with  stochastic  resistance  and 
deterministic  stress.  This  means  that  a  structure  will  fail  if  its  resistance  R  is  below  a 
constant  failure  level  s. 

In  time  the  strength  of  a  structure  will  degrade  due  to  deterioration.  Since  deteri¬ 
oration  mostly  is  observed  according  to  a  particular  time-unit,  we  subdivide  the  time- 
axis  (0,  oo)  into  time-units  ((n  -  1)A,  nAj,  with  n  =  1, 2, 3, . . .  and  time-unit  length  A. 
Often  a  structure  is  planned  to  function  for  a  finite  time-interval,  say  (0,NA].  If  the 
structure  has  met  its  needs  it  will  be  replaced  at  IV A.  In  general  this  replacement  will 
be  an  improvement.  With  the  passage  of  time  new  techniques  come  available  and  can 
result  in  a  higher  original  resistance  and  a  lower  deterioration  rate.  In  our  terminology 
an  improvement  differs  from  the  maintenance  action  preventive  repair. 

'  Suppose  that  in  time-unit  n  the  structure’s  resistance  suffers  a  stochastic  degra¬ 
dation  Dn,  D„  >  0,  and  assume  the  resistance  at  time  zero  equals  r,  where  r  >  s. 
Consequently,  the  resistance  in  time-unit  n  can  be  written  as 

R.  =  r-£tt,  »  =  1 . N.  (1) 

l*>  1 


See  also  figure  1. 

We  study  a  physical  deterioration  process,  where  the  expected  resistance  is  linear' 
decreasing  with  time.  Or  in  other  words:  the  expected  ainouuts  of  damage  per  time- 
unit  are  equal  to  each  other.  Examples  within  the  field  of  hydraulic  engineering  are 
according  to  Vrijliug,  Klattcr  &  Kuiper  (15)  settlement1  and  the  behaviour  of  a  dune 
under  the  continuous  action  of  waves  and  wind.  Other  examples  arc  corrosion  of 
uuderwater  pipelines  and  the  scour-depth  near  underwater  footings  of  a  bridge  (Tang 
[personal  communication]). 

Instead  of  supposing  independence  we  make  weaker  assumptions:  i.e.  exchange¬ 
ability  and  l\  symmetry. 

'When  &  civil  structure  is  built  on  clay  some  water  present  in  the  day  will  be  pressed  out  with 
passage  of  time,  causing  sinking  of  the  structure. 


The  vector  Dy  of  N  uncertain  amounts  of  damages, 

Dy  =  (2) 

is  assumed  to  be  exchangeable.  Dy  is  exchangeable  if  the  probability  density  is  invariant 
under  all  JV!  permutations  of  the  coordinates,  i.e.  if 

Pr  {A  =  ,  Dn  =  6n  }  =  Pr  {A  =  6*  1, . . . ,  Du  =  6*y} ,  (3) 

where  tt  €  N\  is  any  permutation  of  l,...,iV.  The  variables  of  an  infinite  sequence 
Doo,  where  D*,  =  limy-.^  Dy ,  are  exchangeable  if  A>  •  •  • ,  A»  are  exchangeable  for 
each  n  (see  Mendel  [10]). 

If  the  coordinates  of  Dy  are  not  exchangeable  Mendel  [9]  provides  a  way  to  trans¬ 
form  them  in  such  a  way  that  they  are  exchangeable.  Also  if  the  expected  accumulated 
degradation  is  not  linear  in  time,  but  exponential  for  example,  it  would  be  possible  to 
apply  exchangeability.  In  a  subsequent  paper  we  will  discuss  also  nonlinear  expected 
deterioration. 

In  all  examples  it  is  assumed  that  the  decision-maker  has  subjective  information 
about  the  average  amount  of  damage  determined  with  regard  to  the  time-horizon 
(0,  JVA].  Consequently  we  condition  the  probability  of  6i , . . . ,  0y  on  the  average  amount 
of  deterioration  over  (0,  JVA],  i.e.  on  EtLi  A  ~  Given  a  constant  value  of  the  sum 
the  probability  density  of  the  underlying  amounts  of  damages  is  furthermore  assumed 
to  be  uniform  or 


Pr  |  A  »  . . . ,  Dft  -  6ff 


for  (6»,...,6y),  €  11+ .  A> ...»  Dy  arc  now  said  to  be  symmetric2  or  fi 

isotropic  (see  Mendel  [10]  for  details).  This  leads  to  a  distribution  uniform  on  simplices. 

The  likelihood  of  realisations  6\ , . . . ,  6n>  1  <  n  <  W,  given  that  Dy  is  k  symmetric3, 
can  be  obtained  by  integrating  the  uniform  distribution  out  on  the  simplex 

EA-M  (5) 

ta| 

2tp  symmetric  stochastic  variables  ....  Ds  have  a  uniform  distribution  on  Df  =»  N&  (see 
Mendel  (lOj). 

3li  symmetry’  imjUica  exchangeability,  on  the  other  hand,  an  exchangeable  measure  which  is  I) 
symmetric  in  two  coordinates  is  h  symmetric. 


=  Pr  |a  —  ,Dy  » 


G 


over  the  (n  +  l)th  through  the  Nth  amount  of  damage.  This  can  be  achieved  by 
application  of  the  Dirichlet  integral  and  results  in 


(iV-l)n 

(N6)» 


Sift' 

N9 


N-n-l 

/[0, *«](£?=!  6i)t 


(6) 


where  (JV  —  l)n  is  the  usual  shorthand  for  (N  - 1)  •  •  •  (N  - 1  -  n  + 1).  /[„, 5](x)  =  1  for 
x  €  [a,  b ]  and  equals  to  zero  elsewhere.  For  a  proof  of  (6)  see  Mendel  [10]. 

An  interesting  case  arises  when  taking  the  limit  for  N  to  infinity.  Then  we  obtain 
a  product  of  n  exponentials, 


=  (7) 

This  result  coincides  with  the  classical  Bayesian  exponential  model,  where  the  coordi¬ 
nates  of  Dn  arc  mutually  independent  and  identically  exponentially  distributed  random 
variables  with  mean  9. 

The  probability  density  of  ^,...,6*  can  be  found  by  formulating  a  finite  de  Finetti- 
typc  representation 

P(b i*  •  •  • » tf*)  -  r°  *  •  •  • » (8) 

/piiO 

where  ir(9)  is  a  prior  density,  which  is  uniquely  determined  by  tho  probability  density 
of  6i, . . .  ,Sn  (see  de  Finetti  [7]  and  Barlow  &  Mendel  (3)).  Mendel  [10]  provides  a  way 
to  obtain  a  prior  for  the  infinite  case.  He  gets  the  inverted  gamma  distribution 

*m  =  —j»*<r+i,cxp{-~]  w«),  0) 

fi,r  >  0,  where  the  parameters  fi  and  r  may  be  based  on  data  and/or  subjective 
opinions  of  experts.  Au  advantage  of  this  prior  is  that  the  calculations  lead  to  closed- 
form  results  iu  both  the  finite  and  infinite  case. 

By  using  the  prior  (9)  in  the  representation  (8)  we  get,  after  applying  the  binomial 
formula,  for  finite  N 


. «.)  = 

_  (W-l).  *Z?'(  N-n-l 

N‘  h  V  * 

=  (Jj  -  1).  "y  ‘  (  A1  -  »  - 1 

N"  a  l  ■ 


r(*  +  n  +  r) 
M**r(r) 


(10) 
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and  for  taking  the  limit  of  N  to  infinity 


JimM, 


r(n  +  r)  £ 

’  no  M"+r 


(ii) 


With  (6),  (9)  and  (10),  for  finite  iV,  and  with  (7),  (9)  and  (11),  for  infinite  N,  we 
can  derive  the  posterior  density  of  the  expected  amount  of  damage  per  time-unit  by 
updating  our  prior  with  available  data  after  applying  Bayes’  theorem: 


J)(5l, . .  • ,  ^m) 


(12) 


Note  that  for  the  infinite  case  the  posterior  is  again  an  inverted  gamma  distribution 
with  parameters  p  +  Ef=i  and  r  +  m.  That  is  why  the  inverted  gamma  distribution 
is  said  to  be  natural  conjugate  to  our  likelihood.  For  calculations  in  detail  see  Cooke, 
Misiewicz  &  Mendel  [5j.  They  also  deal  with  censored  data.  We  get  this  kind  of  data 
in  the  case  of  imperfect  inspection,  i.e.  the  actual  resistance  will  be  determined  with 
uncertainty. 
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3  Decision  problems  from  hydraulic  engineering 

In  this  section  we  will  apply  the  theory  developed  in  the  former  section  to  two  examples 
taken  from  the  field  of  hydraulic  engineering.  Subsequently  the  decision  variables  per 
act  will  be:  the  preventive  repair  interval  with  no  inspection  (§  3.1)  and  the  inspection 
rate  with  given  preventive  repair  level  (§  3.2). 


3.1  Preventive  repair  interval  with  no  inspection 


In  some  cases  a  decision  maker  needs  to  know  when  he  has  to  carry  out  a  preventive 
repair,  while  there  is  no  need  or  no  possibility  for  inspecting  the  structure.  The  ex¬ 
pected  preventive  repair  time  is  fully  determined  by  a  certain  safety  norm,  laid  down 
by  the  government.  An  example  is  the  minimal  permitted  probability  of  inundation 
of  a  region  protected  by  a  dyke  ring  in  the  Netherlands  (see  [6J).  Assuming  that  the 
cost  of  preventive  repair  is  constant  in  the  sense  that  it  is  fully  determined  by  having 
the  maintenance  people  and  the  right  materials  at  your  disposal  and  uot  by  the  actual 
amount  of  degradation  (provided  that  it  has  not  been  failed),  so  a  dyke  manager  wants 
to  postpone  preventive  maintenance  as  long  as  possible. 

Suppose  Pnofw(n)  is  the  safety  probability  norm  provided  for  each  time-interval 
(0,  nA],  1  <  n  <  N,  where  />«,**,(»)  is  strictly  aondecreasiug  in  n.  The  probability  of 
failure  should  not  exceed  P«^(n)  or  with  (1) 

Pr{ft.<s}  =  Pr|fjDl>r-„|<ft1».(n),  ,,  =  l . N.  (13) 


Hence  a  preventive  repair  should  bs  carried  out  just  before  that  n,  where  the  probability 
of  failure  exceeds  the  safely  norm. 

The  probability  of  failure  in  time-interval  n  can  be  calculated  with  (7)  as  follows: 


with  r  -  s  <  NO  and  * (0)  the  inverted  gamma  distribution. 

With  theorem  3  of  the  Appendix  it  follows  that,  for  the  finite  case, 
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(15) 


-'-fifT-.'H'- 


r  —  s]w-<  |Y  —  sl,_1 

~W. 


n  =  1,...,JV,. 


-1 

l  m  1  ’ 

. ,  r  -  s  <  NO,  and  with  theorem  8  of  the  Appendix 
lim  Pr  (y'A  <  r-  5  -^7  V  A  =  #1  = 

(-V!- 


n  1 

=  i -Y — z — 

k(i- D! 


r  -  s 


0  J 


i-i 


exp 


(16) 


n  =  1,2,...,  for  N  -*  oo. 

(15)  is  a  beta  distribution  in  (r  —  s)/N9  with  parameters  n  and  N ~n  (see  theorem 
3  of  the  Appendix).  The  latter  is  the  gamma  distribution  in  (r  -  s)/0  with  parameters 
n  and  1,  which  is  exactly  the  distribution  of  the  sum  of  n  mutually  independent  and 
identically  exponentially  distributed  variables  with  mean  0 . 

The  mean  of  A,  the  amount  of  damage  in  time-interval  (Q,nA]»  equals  nO 
in  both  the  finite  and  infinite  case.  As  we  sec  the  expected  resistance  is  indeed  linear 
decreasing  in  time.  The  variance  of  A  in  case  of  the  beta  distribution  (finite 
time-horizon)  is  n02,  while  the  variance  equals  nO2  for  the  gamma  distribution 
(infinite  time-horizon).  The  variance  for  finite  N  is  zero  if  n  -  N,  which  is  obviously 
correct  since  we  are  conditioning  on  jf  A  -  0 . 

For  finite  N  the  discrete  lifetime  probability  density  function,  given  0,  is  a  binomial 
distribution  (see  theorem  2  of  the  Appendix), 


Pr{x>Szn£>>*  =  = 

t  M  I»1  'v  l*\  J 


for  n  =  1, . . . ,  JV,  r  -  s  <  NO  and  with  parameters  (r  -  s)/N0  and  Ar.  For  iufiuitc 
N  this  lifetime  probability  density  approaches  a  Poisson  distribution  (see  theorem  0  of 
the  Appendix): 


Jim  P,  {  g  D,  S  r  -  ,  "  g  D, ->  r  -  S  |  i  g  D,  =  «}  = 

1  fr-sl"-*  f  r  -  si 

“ornjirn  ■*{-.- H- 


(IS) 


for  n  =  1,2,...  and  parameter  (r  -  s)/0.  The  mean  life  is  1  +  -~(r  -  s)J9  and 
1  +  (r  —  s)/9  respectively. 
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R 


0  A  2A  3A  4A  5A  6A  7A  8A  9A  10A  •••  NA 


t  t  t 

Inspection  Inspection  Inspection 


Figure  1:  The  (discretised)  resistance  Rn  in  time-unit  n  with  inspection  interval  length 
k  =  3.  At  the  first  inspection  (t=3A)  i?3  >  p,  hence  no  preventive  repair  takes  place. 
At  the  second  inspection  s  <  R$  <  p,  so  there  is  no  failure,  but  a  preventive  repair  is 
needed.  The  total  cost  is  2c/  +  cp.  Suppose  k  =  9  then  there  would  have  been  a  failure 
in  the  first  inspection  interval  with  cost  cp. 
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3.2  Inspection  rate  with  given  preventive  repair  level 

In  structural  engineering  corrective  repairs  and  the  collateral  losses  are  too  expensive 
and  consequently  need  to  be  avoided  by  well  planned  inspections  and  preventive  repairs. 
In  this  example  inspection  is  needed  to  know  whether  we  need  to  carry  out  a  preventive 
repair.  Object  of  study  is  to  obtain  that  inspection  rate  leading  to  minimal  expected 
cost.  A  high  inspection  rate  will  avoid  a  failure  (almost  surely),  but  entails  high 
inspection  cost.  No  inspection  at  all  leads  to  failure  cost,  so  the  optimum  inspection 
rate  lies  in  between  these  two  extremes. 

Inspection  will  be  executed  with  a  constant  rate:  with  an  inspection  interval  of  k 
times  A,  where  k  can  be  chosen  from  1, . . . ,  N.  Hence  the  inspection  times  form  the 
set 

jkA,  j  =  1(1)*,  (19) 

where  K  —  [N/k\  is  the  greatest  integer  less  than  or  equal  to  N/k.  Furthermore, 
we  are  dealing  with  perfect  inspection ,  i.e.  the  actual  resistance  R  can  be  determined 
without  any  uncertainty.  Inspection  takes  negligible  time  and  does  not  degrade  the 
structure.  Each  inspection  costs  cj. 

During  an  inspection  we  now  have  the  possibility  to  execute  a  preventive  repair.  In 
case  of  no  failure  there  are  two  possible  actions  depending  on  the  inspected  condition: 
no  repair  or  preventive  repair.  We  will  assume  that  a  preventive  repair  takes  place  if 
the  structure’s  resistance  is  below  a  certain  preventive  repair  level,  denoted  by  p,  where 
s  <  p  <  r,  given  bv  the  decision  maker.  A  preventive  repair  will  be  executed  at  the 
end  of  the  jth  inspection  interval,  i.e.  at  time  jkA,  if 

Rjk  ^  P*  (20) 

There  will  be  no  action  if  is  greater  than  or  equal  to  p.  See  also  figure  1.  Future 
research  will  determine  the  preventive  repair  level  p  for  a  finite  and  infinite  time- 
horizon. 

If  the  preventive  repair  level  is  readied  the  structure  need  not  yet  be  in  the  failed 
condition.  However  with  a  low  inspection  rate  a  maintenance  action  can  be  too  late: 
a  disastrous  failure  can  happen.  A  failure  will  be  noted  immediately,  e.g.  inundation, 
while  exceedence  of  the  preventive  repair  level  can  only  be  known  through  inspection. 

Both  the  failure  cost  cp  and  preventive  repair  cost  C/>  are  constant,  where  cp  <  cp. 
Note  that  if  cp  would  have  been  greater  than  or  equal  to  cp  then  the  optimal  value  of 
k  would  have  been  N.  The  cost  function  to  minimise  the  inspection,  preventive  repair 
and  failure  cost  is  given  in  (21). 
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where 


(21) 


c(k )  =  £[jc/  +  cp]pr{#o-i)k>pns  <Rjk  <p}  + 

3= 1 
K 

+  E  t?c/  +  cf]  Pr  {i2o-i)fc  >  p  n  Rjk  <s}  + 
i= i 

+  [Kc/  +  cp]  Pr  {Ria  >  pC\Rn  <  s} , 


Pr  {Rxk  >pC\RN<s}  = 

N 

=  E  Pr  {-Rfffc  >  p  n  -R«-i  >  s  n  i?„  <  s} .  (22) 

n=/0:+l 

Note  that  the  model  is  an  extension  of  the  basic  model  of  Barlow,  Hunter  & 
Proschan  [2]  and  is  also  studied  by  Gijsbers  [8).  The  latter  assumes  a  normal  dis¬ 
tribution  for  R(t)  (in  case  of  a  continuous  model)  and  does  not  give  analytical  results. 
Under  these  assumptions  the  civil  structure  actually  has  a  probability  to  improve! 

Now  we  are  concerned  with  calculating  the  probabilities  in  (21)  and  (22).  For 
tedinical  details  we  refer  to  the  Appendix. 

By  theorem  7  of  the  Appendix  it  follows  that  the  probability  of  preventive  repair 
at  the  jth  inspection,  given  the  structure  has  not  been  preventively  repaired  at  the 
former  inspection  (i.e.  at  time  (j  -  1)A),  conditional  on  0,  equals 


whore  r  -  $  <  NO . 

In  a  similar  manner  we  can  get  with  theorem  5  of  the  Appendix  the  probability  of 
failure  during  the  jtb  inspection  interval,  given  0,  while  a  preventive  repair  was  not 
needed  at  the  (j  -  l)th  inspection: 

Pr  {  flu-.)*  >  P n  Rit  <  i  |  i  £  D,  =  fl)  = 
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r-s  <  N6.  The  probability  of  failure  in  the  time-interval  from  the  last  inspection 
until  the  end  of  the  time-horizon,  given  in  (22),  can  be  calculated  by  summing  the 
failure  probabilities  over  the  time-units  (Kk  + 1) A, . . . ,  JVA , 

As  in  (14)  we  can  take  the  inverted  gamma  distribution  as  a  prior  7r(0)  to  calculate 
the  posterior  probabilities  in  (21).  As  a  result  we  have  closed-form  expressions. 

For  the  infinite  case  the  probability  of  preventive  repair  in  the  jth  inspection  in¬ 
terval,  conditional  on  0,  is  with  theorem  12  of  the  Appendix 


;&Pr  { >  pn  s  <  %  <  p|  ^ ]££>,  =  == 

1  f  fr-s]'*-1  _  (  jk  -  1  \  [P  -  rr  -  pi*-1) 

(j*-l)!\l  9  J  ti  \  *'~1  /  l  9  J  l  0  J  j 


and  the  probability  of  failure  in  the  same  interval,  given  0,  is 


according  to  theorem  11  of  the  Appendix. 

The  optimal  value  of  C(k)  can  be  found  by  minimising  over  k  —  1, . . . ,  JV.  The 
restriction  that  k  >  1  does  not  lead  to  big  problems.  If  the  optimum  inspection  interval 
length  will  be  lower  than  A,  then  our  model  is  wrong  and  we  need  to  subdivide  the 
time-axis  iuto  smaller  time-units. 

As  a  result  of  symmetry  in  the  detrioration  process  we  are  able  to  express  the 
probabilities  of  preventive  repair  and  failure  in  closed-form  results.  Hence  C(k)  can  be 
evaluated  easily  for  every  k. 


4  Conclusions 


In  this  report  we  used  symmetric  properties  of  the  deterioration  process  to  build  main¬ 
tenance  optimisation  models.  A  consequence  of  the  symmetry  is  that  the  expected 
accumulated  damage  to  a  hydraulic  structure  is  linear  in  time.  The  Bayesian  approach 
uses  a  prior  for  the  average  amount  of  damage  with  regard  to  a  finite  or  infinite  time- 
horizon.  The  amounts  of  damage  axe  assumed  to  be  h  symmetric. 

Advantages  of  the  approach  are  that  we  base  our  model  on  an  observable  quantity. 
Instead  of  a  lifetime  or  a  transition  probability  matrix  we  use  the  structure’s  resistance 
suffering  stochastic  deterioration  in  every  time-unit.  Furthermore  it  is  possible  to  use 
data  to  update  our  prior  of  the  average  amount  of  damage,  both  in  case  of  perfect 
and  imperfect  inspection.  We  used  symmetry  in  the  physics  to  simplify  our  model 
considerably,  while  not  assuming  mutually  independent  damages. 

We  applied  the  Bayesian  approach  to  two  decision  problems  from  the  field  of  hy¬ 
draulic  engineering.  The  first  one  is  concerned  with  determining  a  preventive  repair 
interval  given  a  certain  safety  norm.  The  second  is  an  inspection  problem.  How  of¬ 
ten  do  we  have  to  inspect  a  hydraulic  structure  to  ensure  that  the  sum  of  expected 
inspection,  repair  and  failure  costs  is  minimal? 

The  main  advantage  of  the  foregoing  approach  is  that  the  model  is  not  built  on  ad 
hoc  assumptions,  while  the  probabilities  of  preventive  repair  and  failure  in  a  certain 
(inspection)  interval  still  can  be  expressed  in  closed-form  results. 
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A  Technical  proofs 


Theorem  1 


fX  fX-Si  rx- 

/(«,*)  =  [  / 

JSi—  0  */<$2=0  J&m— 1 


*  *  *  /  l<Wn  *  *  *  d6%d6i  =  — r, 

=0  j£.=0  «! 


n  =  1,2,...,  x  €  K+.  (27) 

Proof: 

/(n,  a;)  is  the  classical  multi-dimensional  Dirichlet  integral.  □ 

Theorem  2 

{n-l  n  I  i  ^ 

EA  <sn£D,>z  -£fl,  =  <l  = 

<**> 

n  =5  1, , , . ,  JV,  x  <  NO ,  *  €  Rf,  u/Wch  is  the  binomial  density  with  parameters  x/NO 
and  N,  with  expectation  ^•(x/O)  + 1. 

Proof: 

With  Iff  (Si  ,  62,  •  •  •  1 4|0)  given  in  (6)  it  follows  that 

{w-S  n  1  N  \ 

£A<*n£u(>*  -£a=4= 

/o I  leal  lot  J 

«  /•*-<,  rN*-I2Z 1  ft 

-  LL 

4=o  k'V'-'E&'t,  NO  NO  l  JV0  “  ®°l 

/*  JV  - 1  JV-n  +  lr,  as  1*-*  ,, 


/*  «<  JV  -  1  JV  -  n  +  1  r  as  V 

Jh*o”tJim„l*9  NO  NO  r  jwJ 

JV0  JV0  l1  JV#]  (n-l)! 

/  JV  —  1  \  f  [JCT-1 

V  »  - 1  )  [  “  JVfl]  IJV^J  ’ 


d6n-i  •  •  •  dSi 


where  we  used  the  result  (27)  in  the  last  step  but  one. 


Theorem  3 


P.{t»s.|sE«-<)- 

■  H-AT  I* 


(29) 


n  =  1, . . . ,  N,  x  <  NO ,  x  €  Rf ,  which  is  the  beta  distribution  in  x/NO  with  parameters 
n  and  N  —  n. 

Proof: 

The  proof  follows  directly  from  theorem  2: 


l—l  V  m=l  m=l 


>  X 


1  N  \ 

jy  EC~=# 

iV  ro=l  J 


X 


1-1 


Ml 

n  =  1, . . . ,  N,  x  <  NO,  x  €  R+. 

We  get  the  probability  density  function  by  differentiating  (29)  with  respect  to  x: 


1  fl  X  f  X  1 

B(n,W-n)r  TO  [no\ 


x  1 n_1  1 
NO 


which  is  the  beta  density  with  parameters  n  and  N  —  n. 


P 


Theorem  4 

^0>»*,y) 


rv- sl,‘>  rv-Ei:'* 


r  r 

Ji l«sO  Jif  C«0  "Si*- 1*=0  Js, 

1 


/jiBO 


ld6n  •  •  •  d6i 


n! 


"-e(  <"!)(» -*r-i+,*i-1 


j  <n~  1,  j,n  -  1,2,...,  x  <  y  <  NO,  x,y  e  R+. 


(30) 
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Proof: 

The  integral  can  be  determined  by  successively  integrating  out  the  variables  6n,  <5n_i, 

. . 61. 

J(j,  n,x,y)  = 

J6i=Q  J6j= 0  Jsj+ 1=0  J6n=  0 


__  rs  rx-'Lit=isi  ru-HUi6*  /'«-£*»' '*s>  1 

Jsi—0  4=0  Jsj+i=Q  Js%— j=o  2! 


n-2 


y-X> 

1= 1 


d<5n_2  •  *  *  af^i 


r*  /•*-££ i  s>  l 

ijj=0  JSjsd 

=  -r  -r^ 

JS\—0 


«i= o  (n-i)! 

v-i. 


nn-j 


y-E<5< 


/=! 


1 


-EC* _ 


(y-x]n"i+1d<5i_1...dii  + 


r  ...r* 

•/$i=0  |=i 


-EG*  l 


o  (n  -  j  + 1)! 


;-i  I  n~-»+1 


y-E^ 

1=1 


d<5j_i  •  *  • 


The  first  integral  can  be  solved  by  using  (27)  and  becomes 

(y  ~  z)n“*’fla^“1 

”(n-;+i)!(j-i)r 

By  successively  solving  the  second  integral  we  get: 

(y  -  »)"~i+V~2  (i/-*)""1*1  r*  1 

~  (n  -  j  +  2)!(j  -2)!  (n  -  1)11!  +  4=o  (n  - 1)! 

where  we  used  (27)  (J  -  2)  times.  The  last  integral  leads  to 

(y-*)n  ,  vn 


(31) 


[y  (32) 


r 

J(XZ 


1 


$,=o  (n  - 1)! 


n! 


»!* 


(33) 


By  summing  (31),  (32)  and  (33)  and  multiplying  them  with  n!/»!  we  see  immediately 
that  J(j,n,a:,y)  equals 


1 

u! 


»"-§(")(»-  *rv  =  jj »"  -  E  ( !  j )  or -  *r'+v-‘ 


j  <  n  -  1,  j,w  =s  1,2,...,  a?  <  y  <  JV0,  x,y  e  11+. 


□ 
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Theorem  5 


Pr|£^<  <yfi£4  >y 

[  1=1  i=i  i=i 


*£«-«}■ 

-  (::;)i-a'“{fer-i(::,‘)erfer)' 


j  <11-1,  j,n  =  l,...,N,  x<y<  NO,  x,y  €  R+. 


(34) 


Proof: 

In  evaluating  (34)  we  use  the  likelihood  l  tf  {61,62,  •  •  • ,  6n|0)  derived  in  (6).  The  integral 
can  be  determined  by  first  integrating  out  the  variable  6n. 

Instead  of  the  notation  (N  - 1)  •  •  •  (N  - 1  -  n  + 1)  we  will  use  the  usual  shorthand 
(N  -  l)n.  In  the  last  step  but  one  we  use  J(j,  n  —  l,x,y),  following  theorem  4. 


Pr(x>.< 

(i=l 

=  r  -rs;,,'7 

4=0  J6j= 0  4 

_  r*  rx~T>iml  r 
4t=0  4=0  4 


»£“•*}  - 


<zn£A<yn£A>y 

i=i  i=i 

t. 

fy+i=0 

r"*-Z w*  (N  - 1) 
’  L=v~rr} 


j6*aV-UmX  *1 


6j+i=0 


1- 


EkA 


JW 


at— »— 1 


•  •  dtii 


w-est,**  w 

_  r . .  r°-*  ■ ■  •  r^’" ii -  -Lr-v. . 

y«,=o  V o  ■/«,„=»  (W)«-‘  l  JWJ  *  1  1 

-  (A'-^ [> - &P (— 5i [r' -£(::!) (v-rv-j 


(JW)" 


(  N - 1  \  rt  jLr-JfJLr1  n"1)  nLZLEr"* f JLl ,“1  \ 

[n-ljl'-ml  \lwj  ~£i(i-lJlN0l  IJV0J  J’ 


i  <  n  -  I,  j,n  « 1,. . . ,  N,  x  <  y  <  NO,  x,y  6  II+* 


□ 
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Theorem  6 


j  <  n  —  1,  n  =  1, . . . ,  iV,  x  <y  <  NO,  x,y€H+. 


Proof: 

The  proof  follows  from  theorems  2  and  5: 


i 


<y 

^  i=i 

=  i-S>{£a.  <a?n  >ar 


/~1  vmel 


»i°— I 


+  £  p4EAn<*n£Dm<yn£0w>i,  i  £  A 

<=i+l  1  m=l  mcl  mol  mol 


8 

jw  J  live 


j<n-lJtn-h...,N,x<y<Ndt  xty£  R+. 


Theorem  7 


Proof: 

In  evaluating  the  probability  in  (36)  we  use  the  likelihood  1^(6 1, 5a,  •  •  •  *  <5n|0)  derived  in 
(6).  The  integral  can  be  determined  by  first  integrating  out  the  variable  Thereafter 
application  of  theorem  4  in  the  form  of  J{j,n  —  1  ,x,y)  leads  to  the  required  result. 
Again  we  will  use  the  shorthand  (N  —  l)n  to  denote  (N  -  1)  •  •  •  (N  - 1  -  n  + 1). 


Pr|EA<*n"f;A<»n*<i!><!/  ~I>=4  = 

[  /o  1  i=l  1= il  1=1  J 

_  /*  *<  rv-EL,«( 

J(|cO  JSj  B=0  Js 

m  r  ..  / 

iiicfl  JSj  eO  J6' 


...  / 

/(5j+i=0  Jtf, 

rir-Ef., *  (AT  -  i)n  r  _  Ek.ji, 

i^o  4—sr.>  (no)k  _  jv® 

=  _  [‘  /«-£);> 

Jsi=0  Jij=0  Js  j+i=0 

r-z£‘‘  /v-El.*. 

/fyt=0  JSj+i=Q  A.-,=0  “(OT) 

/  fv  -  i  \  /  r  f 

V»-l/\lAr«J  lAW  J 


W-w-l 


*  *d6i 

liwpr  l  "  iwj  <w-1 ' 
r  r-ltt-  rfl.*  (w - 1)— i  f.  »  T-  ,  ,, 

is,-o  /i,=o  JL-  (N0)n~l  l  -1  1 


/  r,  *  r-  r 
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Theorem  8 


Umpt{|:fl(<x|igA=9}  = 

,  ^  i  r^V-1  f  x\ 

Sr-ousJ  <3ipr#r 

n  =  1, 2, . . . ,  which  is  the  gamma  distribution  in  x  with  parameters  n  and  1/6. 
Proof: 

Follows  from  theorem  3  by  evaluating  the  logarithm  in  the  expression 

[1-j^]W"  =  ®cp{CJV-»)log[l-^]} 
with  a  Taylor  series  and  taking  the  limit  of  N  to  infinity. 


Theorem  9 


>taPr{gA<,nEA>,|iEA  =  «}  = 


i  m1*-1 


(»-l)iW  6i'  ^ 

n  =  1,2,...,  which  is  the  poisson  density  with  parameter  x/6  and  expectation  (x/6)  4- 1 . 
Proof: 

The  proof  follows  from  theorem  2  in  a  similar  way  as  the  proof  of  theorem  8.  Q 
Theorem  10 


j  <  n  - 1,  j,n  =  1,2,...,  x  <y,  x,y€ll*.- 
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Proof: 

The  proof  follows  from  theorem  6  in  a  similar  way  as  the  proof  of  theorem  8.  □ 


Theorem  11 


lim  Pr  A  A  <yn£  A  >y  4]£A  =  4 

j  <  n  -  1,  j,  n  =  1, 2, . . . ,  x  <  y,  x,yeR+. 

Proof: 

The  proof  follows  from  theorem  5  in  a  similar  way  as  the  proof  of  theorem  8. 


(40) 


Theorem  12 


&p4tA<*ni;D(<!/ns<tA<w  j~EA=«} 

*  °°  f  leal  1  Jal  iV  f*l  J 

•  {[I]”' -fi  ( 111"'}  • 

j  <  w  —  1,  jf, n  ss  1, 2, . . . ,  x  <  yf  xty€l4. 


(41) 


Proof: 

The  proof  follows  from  theorem  7  in  a  similar  way  as  the  proof  of  theorem  8.  Q 
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